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We study collective spin oscillations in a spin-1 Bose gas above the Bose-Einstein transition tem- 
perature. Starting from the Heisenberg equation of motion, we derive a kinetic equation describing 
the dynamics of a thermal gas with the spin-1 degree of freedom. Applying the moment method to 
the kinetic equation, we study spin-wave collective modes with dipole symmetry. The dipole modes 
in the spin-1 system are classified into the three types of modes. Frequencies and damping rates of 
the dipole modes are obtained as functions of the peak density. We find that the damping rate is 
characterized by three relaxation times associated with collisions. 



I. INTRODUCTION 

Since the experimental realization of Bose-Einstein condensation (BEC) in dilute atomic gases P, S], there have 
been a large number of theoretical and experimental studies on ultracold atomic gases. Among other things, BEC of 
atoms with internal degrees of freedom attracted much attention. In particular, the dynamics of a spin-1 Bose gas 
has been extensively studied, since the MIT group have succeeded in confining ^■^Na Bose-Einstein condensate in an 
optical trap (so called spin- 1 BEC) H. Such spin-1 BEC exhibits many intriguing phenomena because these spin 
components can exchange each other through spin-spin interaction. Although the strength of the spin-spin interaction 
is typically an order of magnitude smaller than that of the spin-independent interaction, the spin-spin interaction 
have a significant effect on dynamical properties of spin-1 Bose gases [SflSHH 

Most of the studies in the d yna mics of spin-1 Bose gases mainly concentrated on near zero-temperature dynamics in 
a pure condensate [l^, [ll|, [l^.]l3l [l3 . [Tsl . [l6l . . On the other hand, there have also been a growing interest in finite- 
temperature properties of the spinor BEC, taking into account the effect of a thermal component [9l.ll6l. fisl. [l9l. [20l. [2l]| . 
In this paper, we study the spin dynamics of spin-1 Bose gases above the Bose-Einstein transition temperature (Tbec)- 

In fact, spin- 1/2 Bose gases are known to exhibit the collective spin dynamics even well above Tbec |H, H^. In 
JILA experiments a trapped dilute Bose gas of ^^Rb atoms with two hyper-fine states were used to study the 
dynamics of a spin-1/2 Bose gas. Even though the JILA experiments were done in relatively high-temperature regime 
T ^ 2rBECi where the quantum degeneracy has little effect on the thermodynamic properties, this spin-1/2 system 
exhibited collective spin dynamics due to the exchange effect, such as spin segregation [1^, [H, [S^, , and spin- wave 
oscillations There have been many theoretical studies discussing collective spin oscillations in the spin-1/2 



system [26|, |27|, |28|, |29|, |30|, |31l IS^, l33|, |3J1 , which showed very good agreement with the experiments [25|, |35|, [36 

We also expect the spin-1 thermal gas to exhibit collective spin oscillations due to the spin-spin interactiorL|20l. I21j| . 
The purpose of the present paper is to investigate of the collective spin oscillations in a trapped spin-1 Bose gas above 
the Bose-Einstein transition temperature. We extend the work on the spin-1/2 system in Ref. [2^ to the spin-1 system. 
Ref. [2^ discussed the collective spin oscillations with dipole and quadrupole symmetries for the spin-1/2 Bose gas, 
applying the moment method to the kinetic equation [s^, HI] . In this paper we discuss the dipole spin oscillations in 
a spin-1 Bose gas. While in the spin-1/2 system, the components couple each other through the exchange effect, in 
the spin-1 system the components couple through intrinsic spin-spin interaction. We thus expect the spin-1 Bose gas 
to exhibit a richer collective spin dynamics. 

In Sec. II, we first derive a kinetic equation for the spin-1 Bose gas. In Sec. Ill, starting from the kinetic equation, we 
derive moment equations describing spin oscillations with dipole symmetry (29l. [39j. Solving the moment equations, 
we obtain explicit analytical expressions for frequency and damping rate of dipole modes. Asymptotic expressions for 
frequency and damping rate in coUisionless limit and in the hydrodynamic limit are also obtained. 



II. DERIVATION OF THE SPIN-1 KINETIC EQUATION 



In this section, we give a derivation of the spin-1 kinetic equation, by which we will discuss the dipole oscillation 
above Tbec in SeclIIII 

We consider a Bose gas of atoms with spin _F = 1 in an optical trap (spin-1 Bose gas). Each atom has three hyper 
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spin states; mp = 1, 0, —1, so the single-particle state is expressed by a spinor wave function. 



<A(r,i) 



01 (r,t) 
00 (r, t) 



01 (r,t)|l) +0o(r,t)|O) +0_i(r,i)|-l) . 



(1) 



Bold-faced variables are vectors in coordinate space and the single underline indicates a spin variable. The single 
atom Hamiltonian is given by 
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where the double underline indicates a 3 x 3 matrix, and a, (3 indicate spin components x, y, z. The first term describes 
the center-of-mass motion of an atom in a harmonic trap V (r) ~ ^ {'-^x^^ + '^yV^ + '^z'^^)i ^^'^ the mass of a 
single atom. The second term describes the linear Zeeman energy, where B" is the magnetic field of spin-a component, 
and Sf (a = x, y, z) denotes the spin-1 matrix of spin-a component 
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and g is Lande's g factor and fis is Bohr magnetron. In the third term, we have introduced a generalized form of an 
external field B^^^ that is coupled to quadrupole operator Q"^ , which is defined by 



(4) 



In the spin-1 system, quadratic Zeeman effect described by zz-component i?^^ plays an important role, as observed in 
dephasing and quenching phenomena [3,[lB,[l3 the presence of an external bias field. In Eq.Q, we have introduced 
more general form of the tensor field B"^ , which can express the quadratic Zeeman effect in the presence of magnetic 
field in an arbitrary direction or some kind of two-photon coupling field. In Sec. III. we will only consider the effect 
ofB|^ 

Next, we consider the effective interaction for two F = 1 atoms. The two-body interaction is modeled by a (S-function 
pseudo potential, which is described in the form 
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where the coupling constants go and g2 are given by 
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with ao and a2 being the s-wave scattering lengths for collision channels with total spins J- = and 2, respectively. 
In second quantized form, the many-body Hamiltonian of this system is given by 
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where (r) is the Bosc field operator satisfying the commutation relation 
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and i,j indicate the hyperfine spin states nip = 1,0, —1 and the hat (") indicates a second quantized operator. 

We next define the time evolution of the system in terms of the statistical density operator, and introduce the 
Heisenberg representation for the field operators. The dynamics of the system is described by the statistical density 
operator p (t) , from which one can obtain the expectation value of an arbitrary operator O (which has no explicit 
time dependence) 



(0)t =trp(t)0. 

The state of the many-body system evolves in time according to 



^ dt 
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With Eqs.(l9]) and ((TO)) . the equation of motion for the quantity {0)t is given by 

d{d)t 
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It is convenient to introduce the time evolution operator U {t,to), which obeys the equation of motion 
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with U {tQ^to) ~ 1. Here tg is the time at which the initial noncquilibrium density matrix p (to) is specified. One can 
then express the time evolution of p (t) as 



i3{t)=U{t,to)p{to)U^ {t,to). 
Thus, the time evolution of {0)t can be written as 

{d)t ^trU {t, to) p (to) {t, to) 6 = tvp (to) {t, to) 6 U {t, to) = {6 (t)), 



(13) 
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where O (t) = Z^t {t, to) 6 U {t, to) is the operator in the Heisenberg picture, which obeys the Heisenberg equation of 
motion 
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The Heisenberg equation of motion for the Bose field operator is given by 
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A. The general kinetic equation of the nonequilibrium system 
In order to describe time evolution of the nonequilibrium system, we introduce the Wigner operator 

% (p, r) = y" dr'e'P-"-'/'-^*] (j. + j.72) (r - v' /2) , 

and define the semi-classical distribution function 

(p, r, t) = I ^(p, r) I j) = trp (t) W,, (p, r) , 



(17) 
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where p {to) is the statistical density operator. Knowledge of this function allows one to calculate various nonequilib- 
rium physical quantities, such as the local density given by 
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where 



{i I n (r, t) I j) = (r, i) = (^Pj (r, t) (r, <)) . 

From this local density, the total density, magnetization and quadrupole moment are given by 

n{r,t) = Tr ^(r,t) , 
M" (r,t) EETr£"n(r,t), 
A"f^ (r,t) = Tr Q"''n(r,t), 



(20) 

(21) 
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where "Tr" denotes the spin trace. 

In order to derive a kinetic equation for Wij (p,r,t), it is useful to separate the Hamiltonian into a part that 
describes the mean-field dynamics Hmf and the remaining part H' 



H (t) ^ Hmf (t) + H' (t) , 

where the leading mean-field Hamiltonian is 

^MF (t) = E / ^^^I W \Hgir,t) + U^ir,t) I (r) , 
with [/]_ (r, i) being the Hartree-Fock (HF) mean-field potential 

(r, t) = goR (r, t) + g^n (r, t) 1 + 52 ^ S^m (r, t) + 52 ^ (r, t) 

a a 

The perturbation contribution in H' (t) is 
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where for convenience, we have defined the following matrix describing the two-body interaction 

Vu'jj' = gaSij'Si'j + g-iSij ■ S^j,. (28) 
In writing Eq. ([M)l . H' [t] is viewed as a perturbation to the zeroth-order Hamiltonian HyiF {t). Noting that Ui (t) 

is the self-comsistent HF mean-field, Hmf (t) defines excitations of the system at the level of the time-dependent HF 
approximation. It will be shown later Hmf {t) contributes to the free-streeming term of the kinetic equation, while 
H' (t) accounts for binary collisions between atoms. 

With the above definitions, the equation of motion for Wij (p, r, t) is written as 



d_ 
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= -trp {t) \W^, (p, r) , Hmf {t)\ + -^tvp {t) [Wy (p, r) , H' (t) 



(29) 



For calculating the first term of Eq. (|29[) . we assume that the macroscopic variables are slowly varying in space, and 
make use of the approximation 



Uij (r - r', i) « riij (r, t) - r' • Vr^y (r, t) . 

We also use analogous approximations for V, B", and B^^. 

With these approximations, we obtain the kinetic equation for the distribution function 
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where [ , ] and { , } represent the commutator and anticommutator for the 3x3 matrices. The left hand side of 
this equation represents the free-streaming term which is derived fr-om the first term on the right hand side of Ea. (|29p . 
while the right hand side of this equation represents the collision term which is derived from the last term of Eq. ([29|l . 
respectively. We have defined the effective potential using Eg. ipS]) as 



air,t) = Uoir,t) + lhir,t), 



(32) 



where 

Ug (r, t) = V (r, t)l + gfiB J2 ^" ^) + E Q"^ ■ (33) 

The term on the right-hand side of Eq. (j3ip is the collision term, which is given by 

I I I i) - ^TrP (io) \w,j (p, r, t) , H' (t)l . (34) 
- in I J 

The reduction of this term to the form of a binary collision integral is a lengthy exercise. In Appendix [Aj we provide 
a detailed derivation of the collision integral using perturbation technique. The final result is written as 
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is the collision integral and 5Un is the second-order effective potential. Explicit expressions of these 



where 

terms are given in Eqs. (jA18[) and (|A19p in Appendix [XI For a dilute Bose gas considered in this paper, one can 
neglect 5^ [3ll,|4|. 

In summary, we have derived the kinetic equation for the distribution function ^ for the spin-1 Bose gas: 
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where the total effective potential U. is defined in Eqs. ([26|) . ([32]) . and p3|) . and explicit expression for the collision 

integral is given in Eq. (jA20[) . The above kinetic equation is one of the main results in this paper. Starting from this 
kinetic equation, we will discuss the collective modes of the spin-1 system in Section III. 

III. MOMENT METHOD FOR A TRAPPED SPIN-1 GAS 



In this section, we consider small amplitude spin oscillations around the fully polarized state where all the spins 
point up, which implies that all atoms are initially in the state | F = l,mF = 1). We linearize the kinetic equation 
Eq. ([36|) around this initial state. The initial distribution function m? (p, r) is given by 



W^i°i(p,r) = /o(p,r), 
W° (p, r) = otherwise. 



(37) 



Here we assume that the temperature is approximately twice that needed for Bose-Einstein condensation TbeCj and 
thus the initial distribution /o (p, r) in Eq. ([37|) is given by the Maxwell-Boltzmann distribution. Moreover, we assume 
that giu/k^T ^ 1 and thus the mean-field does not affect the center-of-mass motion. With these assumptions, the 
initial distribution in Eq. (j37p is given by 
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1/2 

where Ath = (27r?i /mk^T^ is the thermal de Broghe wave length, f3 = l/k^T is the inverse temperature, //q is 
initial chemical potential, and uq (0) = e^^° /X'^^^ is the initial density at the center of the trap potential. We then 
substitute ^(p,r,t) = ^°(p,r) + i5^(p,r,t) and £(p,r,t) = £° (p, r) + J^(p,r, i) into Eq.® to obtain the 

linearized kinetic equation. Since we assume T ^ TbeC: we neglect the Bose enhancement factor taking —^ Sij 
in the collision integral Eq. (jA20[ ). 

For simplify, we assume that the uniform magnetic field is applied in the z direction. Then, the initial effective 
potential is given by 

£0(r) = y(r)l+eor + ei2'' 

+ Wong (r) + Wono (r) I + W2Y, (^) + W^2 ^'^0 (r) , (40) 



2 2 i— I 

where the linear Zeeman energy is eo = gii^B. The quadratic Zeeman energy ei = -^j[^B^ corresponds to B^^ [3|, 
where vm is the hyperfine splitting. As noted above, we consider the high temperature region (7,710 (0) /fesT ^ 1, 
where the gradient of the effective potential can be approximated as 



Vr£°(r) w VrF(r)l + Vreo£^ + Vrei( 
VrV(r)f. 



(41) 



Since the external potential and magnetic field are assumed to be time-independent, the fluctuation of the effective 
potential is given by 



5^ (r, t) = (r, t) + go5n (r, t) 1 + 52 (r, + 52 ^ 5M" (r, t) 



(42) 



where we have defined the fluctuation of magnetization as (5Af"(r,t) = Tr^"(5n (r, t). Considering the high- 
temperature regime g,;no(0)/fcBr <C 1, and assuming small amplitude fluctuations 5n <^ no, we approximate 
VrSa » 0. 

With these approximations, we arrive at the linearized kinetic equation for 6W . 
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where we have introduced simplified notations for several matrices, which are defined in Appendix [B) The linearized 
collision integral dSW/dt\^^ii is given in Appendix [Cl 

We take moment of the linearized kinetic equation Ea. (|^5)) for general quantity Xa (p, r) [1^, [13, H^- We define the 
expectation value of the fluctuation from the initial state and the part of this which depend on the collision term for 
Xa (P, r) as 
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Using the above moments, we can derive the general moment equation for {xa (p, r)) from the linearized kinetic 
equation Eq.([321)- 



7 



In general, the collision terms couple several components of the matrix (xa), and thus the density n, magnetization 

M" and quadrupole moment A"^ are coupled to each other. We thus obtain form moment equations in matrix. With 
the appropriate chose of the function Xa (p,r), one can discuss several collective modes for the spin-1 Bose gas. In 
this paper, we consider oscillations of the following three quantities: xo = 1; Xi = Xi = Px/ni.- 



A. Precession mode 



First, we consider the precession mode, taking xo = 1- Inserting xo into the general moment equation, we obtain 
the moment equation: 



;^(Xo) - ^eo(Xo)^ - ^^iiXo)_^ - »^mf(Xo)'" = Uo>_™", (47) 



where we have introduced simplified notations for several matrices, which are defined in Appendix [B] The mean-field 
frequency is defined as 



(j^) ^ gj^"-o (0) 
2^2 ' 



- '-^i^. (48) 



where !F denotes the total angular moment or 2. To obtain a closed set of equations, we must truncate the collision 



terms (xo)'^°" by expanding the fiuctuations in the distribution function 5Wij in powers of position and momentum. 



For precession modes and dipole modes (which we will consider in the next subsection), it is sufficient to take the 
following form: 

5W,, (p, r, t) ^ fo (p, r) [aj,^' (t) + (t) x + a'i {t) p,] . (49) 

The coefficients in the expansion can be related back to the set of moments using Eq.([33]); = (xo)ij, ctY — 
f3muj'^{xi)ij, = P{x2)ij- We linearize Ea. (|A20[) in SW (p, r, t) and obtain the collision term as 

/-2(xo)-i-i 3(xo)o-i -6(xo)i-i \ 
iXor" = 270 3(xo)-io 4(xo)-i-i -3(xo)o-i • (50) 
V -6(xo)-ii -3(xo)-io -2(xo)-i-i/ 

Some details on the derivation of the above expression and the definition of the relaxation rate 70 are given in Appendix 

Using the expression for the chemical potential e'^^" ~ N (27r?i) {f3ujho/'2,'rr) with ujho = {i^xOJyLOz expression 
for 7o reduces to 

4 ,,f2.^"' 



7o = -n(0)(a2-ao)^(^— j . (51) 
Here, n (0) is the density at the center of trap given by 

We now solve the closed set of coupled moment equations. Considering the normal mode solutions (x) oc e~'"*, we 



find that Eq. (|47|) is reduced to eigenvalue equations, whose eigenvalues give frequency and damping rate ui ~ VL — i^. 

From Eas. (j47|) and (|50p . we find that the moment equations are divided into three individual sets of equations. 
First set of equations is given by 

t^(Xo)ii ~4i7o(xo)-i-i, 

t^(Xo)oo = 8i7o(xo)-i-i, (53) 
t^(Xo)-i-i = -4i7o(xo)-i-i- 
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These equations describe collective mode involving oscillations of {A^^ mode) . This mode describes purely damped 
motion with uj = —AijQ. This means that N and do not precession. This mode exists even in the absence of the 
magnetic field. The next set of equations in given by 



- ^ - ^) (Xo)io + 2w^p(xo)o-i = 6i7o(xo)o- 
(Xo)o-i = -6i7o(xo)o-i, 



^0 ei o (2) 
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(Xo)-io = -6«7o(xo)-io- 



(54) 
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(M^ ± iMV) and A^^ 
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V2 



These equations describe collective modes involving oscillations of ~ 

(A/^, A^^ mode). Here "+" mode and " mode are degenerate in the damping rate and these frequencies. Precession 
of Af * occurs only in the presence of magnetic field. On the other hand, A^^ modes exist even in the absence of 
magnetic field. The last set of equations is given by 



4uj 



(2) 
MF 

(2)' 
MF 



(Xo)i-i = -12i7o(Xo)i-i, 
(Xo)-ii = -12«7o(xo)-ii- 



(56) 
(57) 



These equations describe collective modes involving oscillations of = {A^^' — A^'^ + iA^^) /2 and A = 

(^j\xx _ j^yy _ iA^y) /2. These two modes are also degenerate in the damping rates and the frequencies have op- 
posite sign. 

We note that the quadrupolc modes exist even in the absence of external magnetic field. Although operators for 
conserved quantities (the total density n and the magnetizations Af") have no values in the precession motion in these 
case, the quadrupolc moments A"^ remain finite, i.e. (xo) ^ since they are non-conserved operators. This result is 
in sharp contrast with the spin- 1/2 system [2^. 



B. Dipole mode 



We next consider collective modes with dipole symmetry in a cigar-sharped trap in the presence of uniform magnetic 
field. To discuss these modes, we need quantities are 

Xi=2^, X2=Px/m. (58) 

Inserting xi ^-nd X2 into the general moment equation, we obtain the moment equations: 

4(xi) - to)^ 

+ ^eo(xi)^ + ^ei(xi)'^ + ^^mf(xi)^ = (xi)""", (59) 
4(X2) + ujI{xi) 



+ ^eo (x2)^ + ^ei (X2)^ + »^mf (X2)' ^ + ^^mf (X2)^ = (X2)' '°", (60) 

where have introduced simplified notations for several matrices, which are defined in Appendix [Bl 

Using the explicit form of SWij in Eq. ((49|) . we obtain the collision terms as (see Appendix ICl for some details) 

'(Xi)ff {xiYif {Xi)Tl{\ /-27o(xi)-i-i 37o(xi)o-i -67o(xi)i-i 

(Xi>8f (Xi>8o' (Xi>8-"i = 37o(xi)-io 47o(xi)-i-i -37o(xi)o-i | , (61) 
.(Xi)iti (Xi>i1[, (Xi>i1'-i/ V -67o(xi)-ii -37o(xi)-io -27o(xi)-i-i, 



u 1 1 ^ Id 

2.S 
2 

i I 

(h 2 4 ID 

FIG. 1: Frequency and damping rate of the dipole modes obtained from Ea. (|69p versus the density at the center of the trap 
potential n(0) for ^^Na. The dashed hues represent the dipole oscillation of the density n whose frequency only depends on the 
trap frequency ujx- The solid lines represent the modes involving oscillations of the superposition of n, M^, and A^^ . 
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where the relaxation rate 70 has been defined in Appendix [Cl and 71 and 72 are defined by 



71 = -n (0) (a2 - ao) (ao + 2a2) 



72 - ^n(O) (ao + 2a2f 



27r 



1/2 



27r 

/3m 



1/2 



(63) 



In summary, we obtained the coupled moment equations associated with xi = a; and X2 = Px/ttl. In addition, the 
collision terms couple different internal components. This implies that dipole modes do not simply describe oscillations 
of the density or magnetization, but oscillations of quadrupole moment are coupled to them. From Eas. ([59)l . ([60]) . (|6ip . 
and (j62p . we find that the dipole mode equations can be divided into five independent sets of equations, describing 
three types of modes. 

First type of modes is involves oscillations of diagonal components of moment matrices, which is described the 
following set of equations: 



' w(xi)ii - i(X2)ii = -2i7o(xi>-i-i, 
w(xi)oo - i(X2)oo = 4i7o(xi>-i-i, 
'^(Xi)-i-i - «(X2)-i-i = -2i7o(xi)-i-i, 

'^(X2)ii + iuj^{xi)ii = ^ (72 + 271 +70) (X2>00 + i 



4i 



- (72 - 271) + 270 



(72 + 271 + 70) (x2)oo + 47o(x2)-i-i 



w(x2;-i-i 



-1-1 = 



(72 - 271) + 670 



2/-1-1- 



(X2> 



-1-1, 



(64) 



Second type of modes involves oscillations of the non-diagonal components, which is described by 

^° '^^-) (Xi)io + w^f(xi)o-i - «(X2)io = 3i7o(xi)o-i, 
(Xi)o-i - «(X2)o-i = -3j7o(xi)o-i, 



h h 

h h 

h n 

h h 



■^MF 
■^MF 



,(2) 



2/10 



1/10 



2uj 



(2) 
MF 



4i I 
— (72 + 271 + 7o) (x2)io + g (871 + 570) (X2)o-l, 

(X2)o-i + «w^(xi)o-i = -7; (472 + 1770) (X2)o-i, 



(65) 



and 



n 

eo (1 



h 



^) (Xi)oi - ^^mf(xi)-io - «(x2)oi = 3i7o(xi)-io, 
(Xi)-io - «(X2)-10 = -3i7o(xi)-io, 



h h 



;2) 
,(0)' 



,(2) 
■^MF 



(X2)oi +«^x(Xl)oi 



, f^O '^l , o (2) 



4i I 

= (72 + 271 +7o) (X2)oi + 2 (871 + 570) (X2)-10, 

(X2)-io + it^x(xi)-io = (472 + 1770) (X2>-10- 



(66) 



Finally, the third type of modes involves oscillations of the off-diagonal component which couple mp = 1 and mp = — 1 
states, which is described by 



,(2) 



Lo — 2— — 2uj-^p 



(Xi)i-i - «(X2)i-i = -6«7o(Xi)i-i, 



,(2) 



(X2)i-i + iujl{xi)i-i = (472 + 871 + 2270) (X2)i-i, 



(67) 



and 



+ 2^ + 2^;^!, 

I 0^0 I (0) I o 



(Xi)-ii - i{X2)-i 



-6i7o(xi)-i 



,(2) 



(X2)-ii +iujl{xi)-i 



(472 -f 871 2270) (X2)-1 



(68) 
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The first set of equations Eg. ([64)) describes dipole oscillations of n, and A^^. Solving the eigenvalue equation 
for the set of (xi)ii, (xi)oo, (Xi>-i-i, (X2>ii, (X2)oo, (X2)-i-i, we obtain the complex dipole mode frequencies 
ti> = 17 — 17 as 



2i 



U! 



y (72 + 271 +70) ± |- 

2i 

— (72 - 271) - 4i7o ± 



1/2 



I (72 + 271 + 70) 



- (72 - 271) + 270 



1/2 



(69) 



In FiglU we plot the frequency (fi) and the damping rate (7) as a function of the density at the center of the trap 
potential. Here we assume ^'^Na atoms, which have antiferromagnetic interaction. We take the following values for the 
physical quantities: iV = 1 x 10^, {Ljj2n,ujy/27r,uj^/2n} = {24, 240, 240} Hz and B 20 x IQ-'^T. For comparison, 
we also calculated the frequency and the damping rate for ^^Rb (ferromagnetic), but we do not plot the results here, 
because distinction between ferromagnetic and antiferromagnetic cases in these modes are almost invisible. 

One can obtain analytical expressions for the frequencies and the damping rates of the dipole mode for Ea. (|69p in 
two limiting cases. For small n (0), or coUisionless limit u; ^ \^\, one has 



(70) 



where the mode frequency only depends on the trap frequency. In general, the magnetic field does not affect these 
modes. For large n (0), or hydrodynamic limit w <C I7I, one has 



±UJx, 



Ai 

-y72, 



UJ : 



0, 



-2«7o, 



(71) 



where we have used the fact that 72 ^ 71 ^ 7o because go w 200(72 in Rb and go ~ 30.g2 in Na. As shown in 
Ea. (|70p . in the coUisionless limit the frequency is degenerate, only depending on the trap frequency. With increasing 
77,(0), the frequencies split and in hydrodynamic limit, reflecting that there are two kinds of modes, i.e., the mode 
involving fluctuations of the density and the mode involving fluctuations of and A^^. The first mode, whose 
frequency only depends on the trap frequency, dose not decay. On the other hand, the second mode involving 
oscillations of and A^^ are over damped. From Eq.fTTj). one can see that in hydrodynamic limit only the density 
oscillations can be excited. 

The second set of equations Eas. (|H5)) and ([5^ describe the dipole mode involving oscillations of M+ and ^4+^ as 
well as those of M~ and A~^. Here these two modes are degenerate in the damping rates while the frequencies have 
opposite sign to conserve the spin density. We only consider the mode for and A~^^. Solving the eigenvalue 
equation for the set of variables (xi)iOi (Xi)o-ij (X2)i0i (X2)o-ii we obtain the dipole mode frequencies uj = ^l — 

as 



1/2 (0) (2) 4i 

^ ^ 2\h^^°^'^^' ^ ^MF + ^MF - y (72 + 271 + 70) 



± 



,(0) 



+ ^MF-f(72 + 27i+7o) 



1/2 ^ ^ o (2) 4* 26i 

2 (^(^o-.i)+3^MF-y72-— 7o± 



,(2) 



Ai 



Aujt 



-MF - y (72 + 270) 



1/2N 



-Atui 



(72) 



1/2N 



In Figl^l we plot the frequency {fl) and damping rate (7) for both ^^Rb and ^'^Na. From Figl^l one clearly sees the 
distinction between ferromagnetic or antiferromagnetic cases. 
In the coUisionless limit, the solutions of Eg. ([72]) reduces to 



OJ : 



(eo + ei) 



Corresponding eigenvectors are given by 

/ (Xi>io 
(Xi)o-i 

(X2)lO 
V(X2)o-l/ 



h 



(eo - ei) ± uJj: 



UJx 


( (Xi>io \ 




(Xi)o-i 








(X2)l0 


i 


V * J 


\(X2)o-l/ 


\ / 



(73) 



(74) 
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FIG. 2: Frequency and damping rate of the dipole mode obtained from Ea. (|72p versus the density at the center of the trap 
potential n(0) for '^^Na (left) and *^^Rb (right). The dashed lines represent the modes having the lowest damping rates among 
the four modes. The solid lines represent the other dipole modes. 



As expected, the quadratic Zeeman energy splits the frequencies of two modes, and the frequencies strongly depend 
on the magnetic field in this coUisionless limit. 



In the hydrodynamic limit, one has 



^ M") I (2) A 4* ^ (2) q- _o (0) f'7K\ 

W W 0, WW [ujIj^^ + j - y72, W W - 3l7o, LO W 2UJI^^ - y72, (75) 

where we made use of the fact 72 ^ 7i ^ 7o, since we now consider ^^Rb and ^"^Na (02 — ao). The last two solutions 
in Eg. ([75]) are strongly affected by the sign of g2- 

In Fig|31 we focus on the mode with the lowest damping rate. The frequency of this mode grows linearly with 
n (0), while the damping rate has a peak. This kind of peak behavior of damping rate is often seen in cross over from 
coUisionless to hydrodynamic modes [2^, [4l| . 

In the coUisionless limit, the damping rate grows as 

2 

Aw « -72, (76) 
while in the hydrodynamic limit, the damping rate falls as 

Aw« (77) 



(472^ +9(4% + w^^)' 



In Fig [21 we also plot these asymptotic behaviors given by Eas. ([7^ and ([77|) . Moreover, we consider how oscillations 
of and are coupled. Calculating the eigenvalue equation to find the eigenvector, we find that this mode 
involves oscillations of both M+ and A~^^ in the coUisionless limit. While in the hydrodynamic limit, the non- 
conserving quantity A'^^ decays quickly, so this mode only involve of oscillations M+. 

The third set of equations Eqs. ([S7|) and describe dipole modes involving oscillations of j4++ as well as A~~, 
which couple mp = 1 state with mp = —1. These modes involve the superposition of quadrupole moments which 
is constructed by multiplying and M^. These two modes are also degenerate in the damping rates and the 
frequencies have opposite sign. We thus only consider the modes for Solving the eigenvalue equation for the 
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Q.2 



UA U.« WS^ 



FIG. 3: The damping rate of the mode with the lowest damping rate in Fig[2] versus the density at the center of the trap 
potential n(0) for ^'^Rb. The solid line is the damping rate obtain from Eq. (|72p and the dashed lines are in the collisionless 
limit Ea. (|76|) and the hydrodynamic limit Ea. (|77p . respectively. 



2m 



3.tM 




FIG. 4: Frequency and damping rate of the dipole mode obtained from Ea. (|78|l versus the density at the center of the trap 
potential 7i(0) for **'^Rb. The dashed lines represent the mode having the lowest damping rate. The solid lines represent the 
other mode. 
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FIG. 5: The damping rate of the mode with the lowest damping rate versus the density at the center of the trap potential n(0) 
at for *''Rb. The solid lime is the damping rate obtained from Ea. (|78|l . and the dashed lines represent the asymptotic behaviors 
in the coUisionless limit and that of the hydrodynamic limit respectively. 



set of (x2)i-i, we obtain the dipolc mode frequencies = — 17 as 



1 / 4 



(0) , r V^J 



,(2) 



4i 



(72 + 271 + 70) + 1270 



± ■ 



'^MF + ^MF - y (72 + 271 + 70) 



,(2) 



1/2n 



(78) 



From Eq. (|75|) . one sees that the oseillations are independent of the quadratic Zeeman energy. In FigHl we plot the 
frequency (fi) and the damping rate (7) given from Eq. (|78p . In the coUisionless limit, the frequency strongly depends 
on the magnetic field as 



(79) 



In the hydrodynamic limit, the mean- field strongly affects the frequency, which can be approximated as 



2uj- 



(2) 
MF 



6«7o, 



y72. 



(80) 



We see that oscillations of the quadrupole moments cannot be easily excited. As shown in Fig|4l one of the two 
modes has a high damping rate, which implies that this mode decays quickly. The both modes have high frequencies, 
which imply that these modes cannot be easily excited. We now focus on the oscillation with the lowest damping 
rate. In Fig|51 we plot the lowest damping rate as well as its asymptotic behaviors. Similarly to Fig|31 we find a peak 
behavior. 

Finally, we briefly discuss experimental possibility of observing spin-wave collective modes predicted in the present 
paper. In the JILA experiment on the spin-wave collective modes of the spin-1/2 system [24| . quadrupole mode is 
excited by quadratically varying frequency splitting A (r) oc , where A (r) is the transition frequency between the 
two states, which is equivalent to the transverse magnetic field in our spin-1 system. In this experiment, frequency 
and damping rate of the spin-wave oscillation was investigated in great detail. In addition to the quadrupole mode, 
the dipole mode of the spin-1/2 system is expected to be excited by a linear inhomogeneous frequency splitting, such 
as A (r) (X 2; [2§|. Analogously to the previous studies on the spin-1/2 system, the dipole mode of the spin-1 system 
can be excited by using linear field gradient. We note that the linear Zeeman field will excite the collective mode 
involving oscillations of the linear magnetization M". In order to excite the modes dominated by the quadrupole 
moment A"^, one needs to apply quadrupole coupling field i?"'^, that couples mp = +1 state and mp ^ —1 state. 
Such an external field may be created by a two-photon coupling field. 
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IV. CONCLUSIONS 



111 this paper, we have studied the coUective osciUation in a dihite noncondensed Bose gas in an optical trap. First, 
we derived the kinetic equation describing the nonequihbrium behavior of the system to discuss the spin-1 Bose gas 
at finite temperatures. Starting from this equation, we derived the moment equation for the dipole mode given in 
Eas. (|59)) and ((60)) . Solving those equations, we found the following three properties of the dipole mode. First, the 
mode can be classified into three independent modes; the first kind of modes involves oscillations for the density n, 
z-component of magnetization and 22-component of quadrupolc moment A^^ , the second kind of modes involve 
oscillations of the superposition of Af+ and and that of and A^^ , and the third kind of modes involve 
oscillations of and that of A~~ . Second, we found that the damping rate is characterized by the three relaxation 
rates. These relaxation rates are expressed as linear functions of the density n (0) at the center of the trap potential. 
Finally, we considered the frequency and damping rate of the dipole mode. It is very interesting that the spin-1 thermal 
gas exhibits strong collective behavior due to the spin-spin interaction even in the coUisionless regime. Moreover, we 
found that A^^ and A^^ modes can exist even in the absence of magnetic field, as shown in Eqs. ([56|) . ([57)) . and (j78p . 
We expect this type of quadrupole modes to be excited by using a two-photon coupling field. 

In this paper, we only considered the normal gas above the BEC transition temperature Tbec- In the Bose- 
condensed phase below Tbec, the spin of condensate and thermal cloud gas will interact strongly For future 

study, we will investigate behaviors of the spin-1 Bose gas including condensate, and study how the condensate 
dynamics is affected by the thermal component. 



APPENDIX A: COLLISION INTEGRALS 

In this appendix, we calculate the collision term Eq. ([34|) explicitly. Following Refs. [Ill, [iOl, we assume that the 
effect of H' [t) is essentially a collision process, which occurs on a time scale much shorter than all other time scales 
in the problem. We also assume that the hydrodynamic variables vary slowly in space and time. With these two key 
assumptions, one therefore expects that the dominant contribution from the commutator of field operators '^i (r, t) 
and terms contained in the effective Hamiltonian involved in Eq. ([M|) to arise from values of r' and t' close to r and 
t. In this situation, we can expand the quantites and J7 that appear in H about this point. It is sufficient to use 



ny (r', t') « (r, t) , £ (r', t') w £ (r, t) 
With this approximation, Hm_f and H' are given as 



MF 



n,,, (r, t) / dv'^, (r') * (r') + n,, (r, t) / dr'*T (r') m, (r') 



(Al) 



(A2) 



^' = ^ E / '^'^^l (^) ^0 (r) 



, (r, t) / dY'% (r') m (r') + n,, (r, t) / dv'% (r') * , (r') 



E/ ^"'jj' 

ii'jj' 

It is now useful to introduce the Fourier transform of the field operator 



(A3) 



(A4) 



where V is the volume of the system. One can then write Hy[-p and H' as 



(A5) 



ii'jj' p 



^ ~ y E E ^"'jj"^Pl+P2-P3+P4'^lpi"I'p2"i'P3'^JP4 



ii'jj' P1P2P3P4 



(A6) 



ii'jj' p 
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It should be understood that the above "local" approximation for the effective Hamiltonian can only be used in 
calculating the collision integral at specific position and time (r,t). 

Using the time-dependent perturbation theory, one can expand a statistical average of an operator O as 

{d)t = tTp{to)\ul{t,to)d{to)Uo{t,to) 



dt'llo{t',to)[ulit,t')d{to)Uo{t,t'),H'{t') Uo{t\to) 



(A7) 



where we have introduced the unperturbed evolution operator 



Uo {t,to) = Texp 



- / dt'HMF {t') 



(A8) 



with " T " being the time-ordering operator. Following Refs. [3l|, |40|, we use the simple approximation for a free 
evolution of a^^ 



ill ^ip^o {t, t') « exp 



-ep {t - t') 



(A9) 



where ep = p^/2m. As usual, we assume that there is no initial correlation so that we can use the Wick's theorem to 
factorizc terms such as 



Within the present assumption of slow variation of the macroscopic quantities in space and time, we can use 

Kpi%P2)to ~ V,P2^y (P'i",i) , 

in the calculation of the collision integral. 

Using these approximations, we obtain the following contribution to the collision term Eq. (j34p as 

trp(io) [w-.j (p, r, , ff' {t) 



j" P1P2P3 *° 
X [E> , (Pi, P2, P3) W<^ (p) ~ E< , (pi, P2, P3) W>^ (p)] 
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VI ^=P3)(*~*')/'' 

i" P1P2P3 

X [E<^. (pi, P2, P3) W>, (p) - S>^. (Pi, P2, P3) W<, (p)] , 



"p+Pi,P2+P3 / e 

i" D,D„D., 



where we have introduced the following notations: 



W< (p) ^ W,, (p) , W> (p) ^ S,, + W.,, (p) , 



(AlO) 



(All) 



(A12) 
(A13) 



S^jj" (Pl,P2,P3) = E ^^rnwV,H"j'r' 
IV rni' i" j' 



W,^^{^?,)W>, {p^)W>„ (P3) 



(Pl>P2.P3) = E ^U'mjV^'^nfj„ 
W rni' j' j" 



+W,'r„MW,'iV2)W> (P3) , 

(Pl)W^^ (P2)W^,MP3) 
+ (Pl)W^/'Z' (P2)W^/i(P3)l , 



(A14) 



(A15) 
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omitted the arguments r and t for simplicity. In the time integral of Eq. (|A12[ ). we set to — oo, which yields 



7t6 (ep + - - ep,) + iV 



Cp ~\~ 6p^ 6p2 ep3 



where T' is the symbol of principal integration. We then obtain 

dW 



I = 



dt 



coll 



mip),sun{p) 



(A16) 



(A17) 



where \coU is the collision integral and 511^ is the second order potential given by 



dm 



SUnip) 



TT 



y ] (5 (ep + Epj ep2 ^pa) ^p+Pi,P2+P3 
P1P2P3 



[{1 +^(P),1< (Pl,P2,P3)} - {^(p),l> (Pi,P2,P3)}] , 

71 E ^ 



P1P2P3 



p -r Cpi £p2 £p:i 



-'P+Pl,P2+P3 



[S^ (Pl,P2>P3) (Pl,P2>P3) 



(A18) 



(A19) 



For a dilute Bose gas considered in this paper, one can neglect 5Un |3l] . |40| . Finally, we replace the momentum 

sum l/VJ2p by the integral J dp/ {2TTTif and the Kronecker delta function VSp^pi by the Dirac delta function 
{2'KKf 6 {p ~ p'). We then obtain the expression for the collision integral 



dm 



coll 



^-^Idp, 



dp2 f dp. 



j(S (ep + ep, - Ep, - Epa) (5 (p + Pi - P2 - P3) 



1 +^(P),1< (Pi,P2,P3)} - {ffi(p),l> (Pi,P2,P3)} 



(A20) 



In 



APPENDIX B: FORMULATION OF THE MOMENT EQUATIONS 

matrices are defined by 



the linearized kinetic equation Eq. (|43p . simplified notations for the 



> 


-5VFio 


-25VK] 







-5 Wo 









1 


-51^10 


-5Wi_ 


mm 



















—577,10 5 


77l_l\ 


-577-01 







-577-11 










-5Wio 


^ 


5Woi 





5Wo-i 


V 


-5W_io 








-5VKio 


-35VKi 


5Woi 





-25 Wo 


35W_ii 


25W_io 






(Bl) 



/ 5/710 + 25770-1 

577,^ = —57701 — 25n_io 
5n-ii 



-5771-1^ 






From Ea. (|43p . one can derive a general moment equation. First, we consider the precession mode by inserting 
Xo = 1 into the general moment equation. The moment equation for the precession mode is given by Ea. (j47|) . where 
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several matrices are defined as 

(xo)io 

Xo)oi -(xo! 

(xo)-io 



(Xo)^ = ( -(Xo)oi (xo)o-i 1 , (Xo)^ = I -(Xo)oi -(xo)o-i 



(Xo^ 






(Xo)io 


2(xo)i 


-(Xo)oi 





(Xo)o- 


2(xo)-ii 


-(Xo)-io 








-(Xo)o-i 


2(xo>i 


(Xo)-io 





(Xo)o- 


-2(xo)-ii 


-(Xo)-io 






(B2) 

/ -(yn)n-i 2(Yn)i-i\ 

D 



We next consider the dipole mode by inserting xi ^ ^ ^-nd X2 = Vxl^n into the general moment equation. The 
moment equations are given in Eqs. (|59[) and (j60p . where several matrices are defined as 



(B3) 



^ (xi)io 2(xi)i-i\ B I ^ ^ 

(Xi) = I -(Xi)oi (xi>o-i , (xi) = -(xi>oi -(xi)o-i 

2(xi)-ii -(Xi)-io y \ (xi)-io 

^ , -(xi)o-i 2(xi>i-i' 

2(xi)-ii -(Xi)-io 



and 



^ (x2)io 2(x2)i-i\ B ( ^ ° 

-2(x2>-ii -(X2)-10 J \ (x2)-io 

, (X2)l0 (X2)l-l\ n / ^ ^^2)10 3(X2)1-1 

(X2)^= I -(X2)oi , {X2f^\ -(X2>oi 2(X2)0-1 

-(X2>-ii / V-3(X2>-ii -2(xi>-io 



APPENDIX C: RELAXATION RATES 



(B4) 



appearing in the linearized kinetic equation 

oil 

Ea. (|43)) . Assuming small amplitude spin oscillations around fully polarized state, we divide W. into the initial part 



In this appendix we consider the linearized collision integral d6W /dt 
a. (j43p . Assuming smal 
and the fluctuation part: 



ffi(p,r,i)=^°(p,r)+5^(p,r,i), (CI) 

where the initial distribution W ° is given by Eqs. ([57)) and ([55)1 . For the linearized form of the collision term, it is 
convenient to express the fluctuation of the distribution function as 

^^(p,r,<) = /o(p,r)^(p,r,t), (C2) 

where fo (p, r) is given in Eq. p8|) . The explicit expression for ij component of Ea. (|C2[) can be related back to the set 
of moments using Eq. ([44|) as 

<Ajj = {xo)v + Pn^^l{Xi)ijX + P{X2)i]Px- (C3) 
Inserting Eg. ljCip into Eq. (|A20|) and linearizing it in 0, we obtain the linearized collision integral dSW/dt : 
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dt ' 



TT 1 



I coll 



^{2TThf,„ 

^ ' LI rnv r J ■] 



^ I '^^^ I / '^^^^ ~ ~ «P3) -^(P + Pi - P2 - P3) 

4/0 (P2) /o (Pa) ^iimjK^^ii 



-2/0 (p) /o (Pi) VJi'yX^in'fJii + FiiH'Vi'iil'Jlj 

+ /o (P2) /o (Ps) VirmjKmj'l'/>i'i' (P2) + yil'm]Vrnilj"4>j"V (P2) 

+ I^mili'^i'lmj'/'i'i' (P2) + Vi7nll'Vii"mj4'l'i" (P2) 
+ V;imjl4iilj"0j"i (Pa) + Vii„yKnij'i(/)j'i (pg) 

+ VirnllVii"rnj4'li" (Pg) + l^im/l l^i'mj0ii' (Pg) 

+ /0 (P) /O (Pl) ["^^i'l, ("^^H"i'i + VU"U')CP^^'' (p) 
(V^'ay + Vivij") (j)j„j (p) 
+ V;;'rnj (Vi'll'l + Viiiui) Sii(l),m' (Pi) 



(C4) 



rates 



Using Eq. ((49|) . one can now express the moment associated with colhsions {xa)if^ in terms of the various relaxation 

(X.)™" ^Y.T. (^o'7o + cf 71 + cf 72) (X/^);™. (C5) 



Im 



In above equation, the sum impHes that the moment equations couple lower moments to higher moments through 

{Xallf''- In addition, the sum X]/m imphes that the collision term mixes different spin states. The coefficients c"^ 
should be determined by taking moments of the linearized collision integral Ea. (|C4p . and the relaxation rates 70, 71, 
and 72 are defined as 



70 



71 



^ / '^^ I '^^ I '^^^ I '^^^ I '^^^ 

<S (ep + Cp, ~ - epj) 6{p + Pi - P2 - P3)/o (p) /o (Pi) 

dr dp dpi / dp2 / dpg 



72 



(27r?i)^ 


Nj 




X 


1 


5152 




TV 




X 


1 




{2nhf 


N J 




X 



x5 (ep 



fpa) (^(P + Pl - P2 - P3)/0 (P) fa (Pl) : 



rfr y dp y rfPi J dp2 J dpg 

x5{ep + ep, -ep2 - Sp.:^) S{p + p^ - P2 - P3)/o (p) /o (Pi) 



(C6) 



(C7) 



(C8) 



We estimate these relaxation rates, assuming the number of atoms = 1 x 10^, trap frequencies 

{ujx /2n,ujy/27r,Lu^/27r} = {24, 240, 240} Hz and the magnetic field B = 20 x lO^'^T. We find 70 = 3.0 x 10"''s~\ 
71 = -6.5 X 10-^s-\ and 72 = 14s-i for ^'^Rb, while 70 = 1.0 x lO^^s-i^ ^ 3 x lO'^s'^, and 72 9.7s^^ for 
^^Na, respectively. 
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